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We theoretically study electronic properties of a graphene sheet on xy plane in a spatially nonuni- 
form magnetic field, B = Bqz in one domain and B — Biz in the other domain, in the quantum Hall 
regime and in the low-energy limit. We find that the magnetic edge states of the Dirac fermions, 
formed along the boundary between the two domains, have features strongly dependent on whether 
Bo is parallel or antiparallel to Bi. In the parallel case, when the Zeeman spin splitting can be 
ignored, the magnetic edge states originating from the n — Landau levels of the two domains have 
dispersionless energy levels, contrary to those from the n ^ levels. Here, n is the graphene Landau- 
level index. They become dispersive as the Zeeman splitting becomes finite or as an electrostatic 
step potential is additionally applied. In the antiparallel case, the n — magnetic edge states split 
into electron-like and hole-like current-carrying states. The energy gap between the electron-like 
and hole-like states can be created by the Zeeman splitting or by the step potential. These features 
are attributed to the fact that the pseudo-spin of the magnetic edge states couples to the direction 
of the magnetic field. We propose an Aharonov-Bohm interferometry setup in a graphene ribbon 
for experimental study of the magnetic edge states. 

PACS numbers: 81.05.Uw, 73.21.-b, 72.15.Gd, 73.63.-b 



I. INTRODUCTION 

Graphene, a two-dimensional (2D) honeycomb lattice 
of carbons, has attracted much attention, because of its 
unusual electronic properties. In the low-energy regime, 
electrons near the two inequivalent valleys, K and K', 
of its electronic structure can be described by mass- 
less Dirac fermionsii^ and exhibit half-integer quantum 
Hall effects,—!^ in contrast to the quantum Hall effect 
of the conventional 2D electrons formed in semiconduc- 
tor heterostructures.'^ The long phase coherence length 
and mean free path, of micrometer order, measured^ in 
high-quality graphenes indicate potential applications of 
a graphene ribbon to coherent nanodevices. 

On the other hand, the electron properties of the 
conventional 2D systems have been investigated in the 
presence of spatially nonuniform magnetic fields. The 
nonuniform fields can cause the formation of charateris- 
tic current-carrying edge states^i^ along the region of field 
gradient, which correspond to the semi VB x B drift mo- 
tion. These states have been referedi^ as magnetic edge 
states in the analogy to the edge states corresponding 
to the E X B drift, along sample boundaries. Their fea- 
tures have been studied experimentally)^ The nonuni- 
form fields can form various magnetic structures such 
as magnetic steps^^ii^ii^ magnetic quantum dotS )-^°'^^i^^ 
magnetic ringsr^ magnetic superlatticeS )^^i^°'^^ etc, and 
play a role of characteristic barriers and resonators of 
electron transport, ^^■^'^ which properties are very differ- 
ent from those formed by electrostatic gate potential. 

A graphene sheet may provide a good experimen- 
tal system for studying the effects of the nonuniform 
magnetic fields as a nonuniform-field configuration can 
be effectively generated by applying a uniform field 
to a curved sheeti^ In addition, the Dirac fermions 
in graphene can have interesting properties under the 



nonuniform fields^^ , which may be different from those of 
the magnetic edge states in the conventional 2D systems, 
since the Dirac fermions have electron-like, zero-mode, 
and hole-like Landau levels, as well as the pseudo-spins 
representing the two sub-lattice sites of the honeycomb 
lattice. Therefore, it may be valuable to study the elec- 
tron properties of graphene in a nonuniform magnetic 
field, which is the aim of the present work. 

In this theoretical work, we study the electronic struc- 
tures of a graphene sheet (on xy plane) in a spatially 
nonuniform magnetic field of step shape, B — Bqz for 
X < and B = Biz for a; > 0, in the integer quan- 
tum Hall regime, based on a nonintcracting-electron ap- 
proach. By solving the Dirac equation, we first inves- 
tigate the low-energy properties of the magnetic edge 
states formed along the boundary [x — 0) between the 
two domains with different fields Bq and Bi, when the 
Zeeman effect is negligible. They are found to strongly 
depend on whether Bq is parallel or antiparallel to Bi. 
In the palallel case of 7 = Si /So > 0, the magnetic edge 
states originating from the n = Landau levels of the 
two domains are dispersionless, contrary to those from 
the 71 ^ levels, where n is the Landau level index. In 
the antiparallel case of 7 < 0, the n = magnetic edge 
states split into electron-like and hole-like levels near the 
boundary. These features, which are absent in the mag- 
netic edge states of the conventional 2D electrons, are 
attributed to the fact that the pseudo-spin of the mag- 
netic edge states couples to the direction of the magnetic 
field. On the other hand, the features of the n ^ mag- 
netic edge states are similar to those of the conventional 
cases. 

We further study the magnetic edge states in the pres- 
ence of an additional electrostatic step potential, V{x) — 
Vb for a; < and V{x) = Vi for x > Q. For 7 > 0, the 
n = magnetic edge states become dispersive, while for 
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7 < an energy gap becomes created around the bipolar 
region in the spectrum of the magnetic edge states. Simi- 
lar features can be found when the Zeeman spin splitting 
is finite, because in the nonuniform field of step shape 
the Zeeman effect behaves as the step potential. 

Finally, we suggest an interferometry setup formed in 
a graphene ribbon for experimental study. In this setup, 
the magnetic edge states can provide partial paths of 
a full Aharonov-Bohm interference loop, therefore the 
properties of the magnetic edge states such as the gap 
of their energy spectra can be investigated. We numeri- 
cally calculate the transmission probability through the 
setup by using the tight-binding method and the Green 
function technique.— i^i^^ The results are consistent with 
the features of the magnetic edge states obtained by solv- 
ing the Dirac equation. We also derive the transmission 
probability, based on the scattering matrix formalism, 
and use it to analyze the numerical results. 

This paper is organized as follows. The magnetic edge 
states are studied without and with the electrostatic step 
potential in Sees. |TT] and lllli respectively. The Zeeman 
spin splitting is considered in Sec. IIVI and the graphene 
interferometry is suggested and investigated in Sec.|Vl In 
Sec. lVn this work is summarized. Throughout this work, 
we ignore the intervalley mixing due to the nonuniform 
field, the validity of which is discussed in Appendix [X] 
In Appendix |B] we derive the transmission probability, 
which may be applicable to other graphene interferome- 
try setups with slight modification. 



II. MAGNETIC EDGE STATES 

We consider a graphene sheet (on xy plane) in the 
nonuniform magnetic field of the step configuration. 



B{x) 



Bqz, X < 
Biz, X > 



(1) 



in the integer quantum Hall regime. Without loss of gen- 
erality, Bq is chosen to be positive, and 7 = Bi/Bq is 
either positive or negative. This configuration may be 

realized with field gradient dB/dx less than 10'* T • A 
and with not-too-strong field strengths Bo and \Bi\ (less 
than, e.g., 100 T). In this case, the mixing between the 
K and K' valleys due to the nonuniform field can be ig- 
nored (see Appendix A), and the electrons in each valley 
can be seperately described by the Dirac equation. 



[vaj ■li + V{x)]^j = E'if, 



(2) 



in the low-energy approximation. Here, J G {K, K'} 
is the valley index, v = •\/3ai/(2fi) ~ 10^ m • s~*, a is 
the lattice constant of graphene, t is the hopping energy 
between two nearest neighbor sites, (Tk = icrx,cry) and 
(Jk' = {—crx,cry) are constructed by the Pauli matrices, 
n = p + eA, p is the momentum measured relative to 
the valley center {K or K' point), A is the vector poten- 
tial, and e(> 0) is the electron charge. The electrostatic 



potential V{x) applied to the sheet must be a smoothly 
varying function of x. The detailed form of V{x) will 
be specified in Sec. IIIII The components of the pseudo- 
spinor j represent the wave functions of the two sub- 
lattice sites (denoted by A and B hereafter) of a unit 
cell of graphene. The Dirac equations for the K and K' 
valleys are connected to each other by a unitary transfor- 
mation U = ffy. This feature leads that the two equations 
have the same energy levels and that their wave functions 
have the relation, ^ k' = U^k- Therefore, it is enough 
to solve the Dirac eqaution for the K valley only. 

Before studying the nonuniform-field cases, we briefiy 
discuss the case of a uniform magnetic field with strength 
B. In this case, the Landau levels are found^ to be 



En = Sgn{n)vy^2\n\heB, 



(3) 



where n — 0, ±1, ±2, • • • is the Landau level index. The 
levels with n > are often refered as electron-like, while 
those with n < as hole-like. In fact, one can obtain the 
Landau levels from the square of the Dirac Hamiltonian, 

El = 2hev^B{m + 1/2 ± 1/2) (4) 

for the K valley, with \n\ — m + 1/2 ± 1/2 and m = 
0, 1, 2, • • • . For later discussion, it is worthwhile to ana- 
lyze E^. The harmonic term with B{m + 1/2) is noth- 
ing but the Landau level of the conventional 2D elec- 
trons, while the next term ±5/2 can be interpreted as 
the effective Zeeman effect of the pseudo-spin; B/2 for 
pseudo-spin up and —B/2 for down. Each Landau level 
with n ^ is composed of the pseudo-spin up (with 
m — |n| — 1) and pseudo-spin down (with m — \n\) states, 
while the n — level comes only from the pseudo-spin 
down (up) states with m = in the K {K') valley. The 
n = level is independent of B as the harmonic and 
effective Zeeman terms exactly cancel each other. 

Turning back to the nonuniform field in Eq. ([1]), we 
consider the case without electrostatic potential, V{x) — 
0, in this section. We choose the vector potential as 
A{x) = Boxy for a; < and A{x) = Bixy for x > 0. 
This choice of A is useful, as the solution of Eq. ([2]) for 
the K valley can be written as '^\.{x) = e~*'^^$^(a;) = 
e~'^''y {(j)j^{x) , (j>*g{x)), where 0j is the wave function of the 
sub-lattice site j € {A, B} and hk is the eigenvalue of Py. 
Hereafter we will measure energy and length in units of 
El (= \/2haPeBo) and I b (= \/ fi/{eBo)) , where I b and Ei 
are the magnetic length and the energy gap between the 
n = and n — 1 Landau levels, respectively, of the bulk 
region with Bq; in these units, the Landau levels are ^yri 
at X < -1 while Sgn{n)^/\j\n at x > l/v/pH- Then the 
equation for 4>j=A,B is found to be 

- + ^) " ^-^]^^-.'«(^) - 0- (5) 

The effective potential^" V;!^ is harmonic. 



\{x + kY + Sj\, x < 
i(7X + kf -f Sj^, X > 0, 



(6) 
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where Sj=A = 1 and Sj=B = — 1- The equation ([5]) has 
the form of the usual Schrodinger equation with poten- 
tial VJfff and eigenvalue j.. Therefore, V^g is useful 

for understanding (j)j^k{x) and E„^k- The form of V^g in 
Eq. ([6]) shows that the pseudo-spin {sj) couples to the 
direction (7) of the magnetic field. 

The solution 0*^^, of Eq. ([5]) for a: < can be expressed 
in terms of the parabolic cylinder functions^^ Di,{z) as 



(a)Y = 2 



(b)y=-1 
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On the other hand, for x > 0, the solution 0^^. is found 
to be dependent on the sign of 7, i.e., on whether Bi is 
either parallel or antiparallcl to Bq, as for 7 > 
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and for 7 < 
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The energy eigenvalue En^k can be obtained under the 
boundary condition of the continuity of the wave func- 
tions in Eqs. ©-(O at x = 0. For the two cases of 7 = 2 
and 7 = — 1, En^k is drawn in Fig. [TJ These two specific 
cases of 7 are enough to understand the characteristics 
of the magnetic edge states. 

The energy levels strongly depend on whether Bi is 
parallel or antiparallcl to Bq. In the parallel case of 
7 > 0, for each n the energy levels gradually change from 
Sgn{n) ^y\n\ to Sgn{n)^/j\n\ as k decreases from positive 
to negative values [see the 7 = 2 case in Fig.[lja)]. This 
feature can be understood from the effective potential 
[see Eq. ^ and Fig. [^c)] . As k increases from nega- 
tive to positive values, the bottom of V^ff{x) moves from 
the region with Bi to that with Bq, passing the boundary 
X = around fc = 0. Therefore, f or k ^ 0, the eigenstates 
have the Landau levels Sgn(n) ■y/7|n| and localize around 
X = —k/"/ (in the unit of Ib), while for fc 3> 0, they 
have Sgn(n)-y/|n[ and localize around x = —k. Around 
k — 0, the two levels of Sgn(n) Y^7|n| and Sgn{n)^/\n\ 
connect smoothly and the corresponding states are lo- 
calized around x = 0. These states have been callediSi 
magnetic edge states and they can carry current along 
the boundary x = 0. 

Contrary to the corresponding states in the conven- 
tional 2D systems, the magnetic edge states in graphene 
have the following different features for 7 > 1. First, 
the edge states with n > behave as electrons, with dis- 
persion dEn/dk < while those with n < as holes 
with dEn/dk > 0. Second, the edge states with n — 
are dispersionless and carry no current, regardless of 
7(> 0). These features come from the nature of the Dirac 
fermions. Especially, the second one can be understood 
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FIG. 1: (Color online) Upper panels: Energy spectra -E„,fe 
for (a) 7 = 2 and (b) 7 = — 1. The energy levels of the n — 
magnetic edge states are highlighted by the filled squares. For 
7 > 0, the n = levels are dispersionless while they split into 
electron-like and hole-like levels for 7 < 0. Lower panels: 
Effective potential V^g'^ix, k) in Eq.® for (c) 7 = 2 and (d) 
7 = — 1. Different values of k = 1.5 (solid line), (dashed), 
-2 (dotted) are chosen. Energy and length are measured in 
units of Ei{= \/2hv^eBo) and Ib{= \/nJ{eBo)), respectively. 



from the fact that for 7 > 0, the effective Zeeman and 
harmonic contributions to E'^^^ cancel each other in the 
both sides of x = 0, as discussed around Eq. (|4]) and 
as shown in the term of V^g representing the coupling of 
pseudo-spin to field direction [see Eq. The case of 

< 7 < 1 can be understood in a similar way. 

Next, we discuss the case of 7 < 0, which is very differ- 
ent from the 7 > case [see Fig.[TJb)]. For large positive 
k, the eigenstates have the Landau levels of Sgn{n)^/\n\ 
or Sgn{n)^y\jn\, while for large negative k (<C 0), their 
energy either increases (showing electron-like behavior) 
or decreases (hole- like). This feature can be understood 
from V/!ff{x,k). As shown in Fig. [Hd), for large posi- 
tive fe, the two local harmonic wells of V^g{x, k) occur far 
from the boundary of x = 0, resulting in the Landau lev- 
els localized in each well. As k decreases to negative val- 
ues, the two local wells move toward x = and become 
to merge into a single well (not harmonic anymore) at 
X — 0, and then the bottom of the merged well increases. 
Therefore, the eigenstates with fc < are magnetic edge 
states localized at a; = and they can carry current along 
X — 0. They are either electron-like {dEn>o/dk < 0) or 
hole-like {dEn<o/dk > 0). One can estimate their group 
velocity for /c <C as {l/h)dEn/dk ~ ztv from the mini- 
mum value of the merged well of V^, where -l-(— ) stands 
for the hole(electron)-like states. In this case of 7 < 0, 
where the magnetic fields Bq and Bi are antiparallcl, the 
eigenstates in n 7^ Landau levels correspond to classi- 
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FIG. 2: The same as in Fig.[T]but in the presence of the step 
potential V{x) with Vo = and Vi = 0.5. For 7 = —1, the 
energy levels are labeled by ul and ur, which are the indices 
of the Landau levels localized at a; ^ —1 and at a; ^ l/\/py[i 
respectively (see text). The n — levels become dispersive 
for 7 > while the energy gap opens between electron-like 
{ur — 0) and hole-like levels {ul ~ 0) for 7 < 0. 

cal motions, so called snake orbits^i^i^ii^, while those in 
the n = level have no clear correspondence to classical 
motions as they have both electron and hole characters. 

For 7 < 0, the effective Zeeman contribution to Ef-^ 
has the opposite sign between the domains of x > and 
a; < 0, as shown in Eq. ([6]). This coupling of pseudo- 
spin to field direction causes an energy barrier at x = 0; 
for example, a pseudo-spin down state has larger effec- 
tive Zeeman contribution at cc > than at a: < 0. As a 
result, as k decreases, the magnetic edge state becomes 
more confined around a; = due to the effective-Zeeman 
barrier. This pseudo-spin feature, which enhances the 
splitting into electron-like and hole-like states, is absent 
in the magnetic edge states of the conventional 2D elec- 
trons. 



III. MAGNETIC EDGE STATES IN AN 
ELECTROSTATIC STEP POTENTIAL 

In this section, we consider an additional electrostatic 
potential V{x) of step shape, 

to the nonuniform magnetic field in Eq. ([1]). Here Vq and 
Vi are constant. This potential gives rise to characteristic 
modification of the n = magnetic edge states, such as 
the creation of energy gap for 7 < 0, as will be seen 
below. Moreover, the modification is directly applicable 
to the case where the Zeeman spin splitting is finite, as 
will be studied in Sec. IIVI 

The step potential V{x) is assumed to be smoothly 
varying in the length scale of the lattice constant a. 
Then, we can still ignore the intervalley mixing and solve 
the Dirac equation in Eq. ^ in the same way as above. 
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FIG. 3: The same as in Fig.g^b) but with (Vo, Vi) = (0, 1.5). 

In Fig. [21 choosing Vq = and Vi = 0.5 (in the units of 
El and Ib), we draw the energy spectra of the magnetic 
edge states for 7 = 2 and 7 = — 1. 

The features of the energy spectra are discussed below. 
For 7 > 0, the eigenstates with large positive (negative) 
k have the Landau levels shifted by Vq (Vi) and do not 
carry current, while the magnetic states around fc = 
have the energy smoothly connecting the Landau levels 
of the large positive and negative fc's. We point out that 
the n — magnetic edge states carry current due to the 
potential step V{x). On the other hand, for 7 < 0, the 
eigenstates with large positive k have the Landau lev- 
els shifted by Vq or Vi, depending on whether they are 
localized in the bulk region of a; < or a; > 0. For conve- 
nience, we introduce the Landau-level indices and nn 
for those localized in a: < and in a; > 0, respectively. In 
this case, the energy levels with n — split, opening the 
energy gap between the electron-like and hole-like mag- 
netic edge states, in contrary to the case without the step 
potential. For 7 = — 1 and Vi — Vq = 0.5, the gap size is 
the same as the step height Vi — Vq. 

We further study the energy gap with varying the 
height V^i — Vo for 7 = — 1. In Fig. [3l we choose 
Vi — Vo = 1.5, which is larger than the energy spacing be- 
tween the n = and n = 1 Landau levels, in contrary to 
the case of Fig. [2Kb) where the height (= 0.5) is smaller. 
In this case, the energy gap of the magnetic edge states 
occurs between the = 1 and nn = —1 Landau levels. 
Moreover, the gap size is no longer the same as the step 
height, but corresponds to the energy difference (= 0.5) 
between the ul = 1 and riR = — 1 levels. 

The above feature can be understood as follows. For 
7 = — 1, the effective Hamiltonian in Eq. ^ is odd under 
the inversion operator 7?,inv, Ti-invX = —x, 

T^inv [v3j ■ n + v{x) ~ -^i±^]7^i„v7^i„v*J 

= - [vdj ■ n + V^(a:) - °^ ' ]7^inv^J 
Vn + V^i 

- (i? - ^li^)7^i„v*J. (11) 

This property is consistent with the facts that the gap 
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(a) Y=2 



(b) Y=2 



(a) 



(b) 



n=0 




(c) Y=-1 

nL=0,nR=0 



(d)y=-1 




FIG. 4: Schematic spin-dependent energy dispersion of the 
n = magnetic edge states in the presence of the Zeeman 
splitting, 2 jV/| — 0.2. We choose the parameters (7, Vo, Vi) 
as (a) (2,0,0), (b) (2,-0.2,0.2), (c) (-1,0,0), (d) (-1,-0.2,0.2). 



center is located at Ec = (Vq -I- Vi)/2 and that the gap 
size is determined by the energy difference between the 
Landau levels just below and above Ec- Similarly to the 
case in Fig. [l] it can be also understood^^ from V^g that 
the states with E > Ec are electron-like while those with 
E < Ec are hole-like. 



IV. ZEEMAN SPLITTING OF THE MAGNETIC 
EDGE STATES 

So far, we have ignored the spin degree of freedom. 
Recently, Abanin, Lee, and Levitov discussed^'^ that in 
graphene the Zeeman splitting can be smaller than the 
Landau energy gap only by the factor of about 10"^, 
due to the exchange interaction, and that it can play an 
important role in the edge-channel transport in the quan- 
tum Hall regime. In this section, we discuss the effect of 
the Zeeman splitting on the magnetic edge states. 

In the nonuniform field B{x) in Eq. ([1]), the Zeeman 
splitting behaves as a spin-dependent step potential. 



V,^{x)=sAzB{x), 



(12) 



where Az = g*jj,B/Ei, s = 1/2 (s = —1/2) for spin- 
up (down) electrons, and /is is the Bohr magneton. The 
exchange enhancement^'^ of the factor can be taken into 
account in g* . We assume that ^ 1 for convenience. 

The Zeeman splitting can be considered as the spin- 
dependent shift of the step potential, V{x) + 
V^{x). The resulting magnetic edge states can be easily 
understand from the features discussed in Sec. IIIII In 
Fig. m we draw schematic energy dispersions of the n — 
magnetic edge states; the extension to n ^ states 
is trivial. For 7 > 0, the n = states are dispersive 
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FIG. 5: Upper panel: Schematic diagram of a graphene- 
ribbon interferometry, which consists of the left (with mag- 
netic field Bo), middle (with —Bo), and right (with Bo) re- 
gions. There is one edge channel along the ribbon edges, 
while there can be magnetic edge states along the left and 
right boundaries, ab and cd, of the middle region, depend- 
ing on the energy of the states and the field configuration of 
{Bo, —Bo). In (a), there appear two magnetic edge states 
along each boundary, and thus the Aharonov-Bohm interfer- 
ence, while no state and no interference exists in (b). Lower 
panel: Aharonov-Bohm effect in the transmission T through 
the interferometry, as a function of the Aharonov-Bohm fiux 
<j!>. The case of (a) is plotted in (c), while that of (b) in (d). 



near a; = 0, as shown in Figs. HJa) and (b). The energy 
difference between the spin-up and down states varies 
from "fAz to Az as k increases. And, the average value 
of the spin-up and down energy levels is shifted by Vq 
{Vi) at large positive (negative) k. On the other hand, 
for 7 < 0, the energy gap of the n = magnetic edge 
states exists even without the electrostatic step potential 
[see Fig.|4l[c)]. For both the cases of 7 > and 7 < 0, the 
sign of the drift velocity (~ dE/dk) is either positive or 
negative, depending on Az, Vq, and Vi. These spin-split 
dispersions show that spin-polarized current can emerge 
between two magnetic domains. 



V. GRAPHENE INTERFEROMETRY 

In this section, we propose an interferometry setup for 
studying magnetic edge states in a graphene ribbon. We 
focus on the case of 7 = —1 in Fig.[2I^b) and demonstrate 
that the energy gap of the magnetic edge states can be 
directly studied by observing the Aharonov-Bohm inter- 
ference of the setup. The interferometry setup is useful 
as well for the other cases of 7. 

We consider a ribbon with armchair edge; a setup with 
zigzag edge will show a similar result. The ribbon con- 
sists of three parts, the current source in the left, the 
middle scattering region, and the drain in the right. A 
nonuniform magnetic field is applied such that B = Bqz 
in the source and drain while B = Bqz ~ —Bqz in the 
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middle region (see Fig. [5]). At the same time, a constant 
electrostatic potential Vq = is applied in the source 
and drain while Vi — 0.5 in the middle region. Then, 
the magnetic edge states form along the left and right 
boundaries ab and cd of the middle region, while along 
each boundary of the ribbon, there is one edge channel, 
which is formed as a mixture of the contribution of the 
K and K' valleysj^i The magnetic edge states are the 
same as those studied in Fig. ^b), when their position 
separates from the ribbon edge by more than the scale 
of magnetic length, Ib = h/ (cBq), so that the overlap 
between the magnetic edge states and the ribbon-edge 
channels is negligible. Note that at each of the scatter- 
ing points a-d in Fig. El^a) , the number of the incoming 
channels is the same as that of the outgoing channels, so 
that the current is conserved. 

The formation of the magnetic edge states depends on 
energy. In the energy range S [0.5, 1] of Fig. [^I^b), where 
each valley supports only one magnetic edge channel, one 
has the edge-channel transport shown in Fig. EJa). In 
this case, there appears an Aharonov-Bohm loop around 
the middle region, which is supported by the two coun- 
terpropagating edge channels along the upper and lower 
ribbon edges_and by the magnetic edge states along the 
boundaries ab and cd. As a function of Bq , one can ob- 
serve the Aharonov-Bohm interference in the transmis- 
sion through the setup. On the other hand, in the energy 
range € [0,0.5], where there is no magnetic edge state 
along ab and cd [see Fig.lSfb)], no Aharonov-Bohm inter- 
ference can be observed. Therefore, the Aharonov-Bohm 
loop can be formed, depending on whether the magnetic 
edge channels exist along the boundaries ab and cd. This 
property allows one to measure the gap of the magnetic 
edge states by varying the energy of the incoming edge 
channel and by modulating Bq. 

We confirm the above proposal numerically by cal- 
culating the transmission probability T{E) through the 
setup, using the tight-binding method^^ and the Green's 
function approach ;^°'^^ Here we skip the details of the 
method and instead refer Refs.—i^. The effect of the 
magnetic field is taken into account by the Peierls phase. 
The strength of Bq_ is set to be about 800 T. At the 
boundaries ab and cd, the magnetic field spatially varies 
linearly from Bq to Bq over the length scales of ?>Ib- We 
choose the ribbon width |a5| and the width |6c| of the 
middle region as about Ibis and 30Zs, respectively, so 
that we can ignore the overlap between the edge chan- 
nels propagating oppositely to each other. 

In Fig. [3 we plot T{E) for E = 0.65 and E = 0.35 as 
a function of the Aharonov-Bohm flux (/> = 2tt\Bq\S/ 4)q^ 
where E is the energy of the incoming states to the setup, 
S is the area of the middle region, and 0o = h/e is the 
flux quantum. As expected, for E = 0.65, one has the 
interference, while not for E — 0.35. This confirms the 
proposal discussed above. 

Finally, we briefly analyze the numerical result for 
E = 0.65 in Fig. UK a). The period of the Aharonov- 
Bohm oscillation is found to be A(j) = 7.8 = 1.24 x 2tt. 



The fact that A4>/{2tt) is larger than 1 indicates that 
the actual Aharonov-Bohm loop has smaller area than 
S, which is reasonable. The lineshape of the interfer- 
ence can be analyzed by using the expression of T in Eq. 
(|B7[) derived in Appendix [B] The lineshape is well fit- 
ted by Eq. (|B7p with parameters oi a — /3 — 0.5 and 
ta = tb = tc = td — 0.23. From this fitting, one can 
get the information of the scattering between the edge 
channels along the ribbon edges and the magnetic edge 
channels along ab and cd. 

VI. SUMMARY 

We have studied the magnetic edge states formed along 
the boundary between the two domains with different 
magnetic fields Bq and Bi in graphene. It turns out 
that the magnetic edge states have very different fea- 
tures from those of the conventional 2D electrons, since 
the formers have pseudo-spin which couples to the direc- 
tion of the magnetic field. As a result, the n = mag- 
netic edge states are dispersionless for 7 = i?i/i?o > 
while they split into electron-like and hole-like current 
carrying states for 7 < 0. The Zeeman spin splitting or 
the additional electrostatic step potential can make the 
n = states dispersive for 7 > and open energy gap in 
the bipolar region for 7 < 0. These features show inter- 
esting manifestation of the Dirac fermions in graphene, 
and the magnetic edge states can play a special role in the 
transport of the Dirac fermions in a nonuniform magnetic 
field, such as spin-polarized current along the boundaries 
of magnetic domains. 

We are supported by Korean Research Foundation 
Grants (KRF-2005-084-C00007,KRF-2006-331-C00118). 

Note added: During the preparation of this 
manuscript, we have been aware of two preprints^i^^ 
where the energy dispersion and current density of the 
snake states in a nonuniform magnetic field of waveguide 
shape are studied. Their results partially overlap with 
our results for the case of 7 < in the Section |TT1 

APPENDIX A: INTER- VALLEY SCATTERING 
IN NONUNIFORM MAGNETIC FIELDS 

In Appendix [XI based on the tight-binding method, we 
show that the mixing between the K and K' valleys due 
to a spatially nonuniform magnetic field can be ignored 
when the field strength and the gradient of the field are 
much smaller than 10^ T and 10^ T • A , respectively, 
which is achieved in current experimental studies. 

We first discuss the matrix elements of the tight- 
binding Hamiltonian of graphene in the presence of an ex- 
ternal magnetic field. For each sub-lattice site j = A, B, 
the Bloch wave functions of an electron is written as 

Rj 



7 



where the sum runs over the potisions Rj of site j, N is 
the number of unit cells, A{ri) is the vector potential, 
k is the momentum of the states, and (j>{r) is the wave 
function of electrons participating in the tt bonding. The 
matrix element ($?|i7|<i>4) of the tight-binding Hamil- 
tonian for the nearest-neighbor hopping is found to be 



,(A2) 



where t is the hopping energy between two nearest neigh- 
bor sites in the tight-binding scheme, the sum runs 
over the nearest-neighbor site pairs of A and B, and C 
denotes the path connecting the site pair. 

Using the matrix element in Eq. (IA2|) . one can esti- 
mate the effect of the magnetic field on the intervalley 
mixing. To do so, we assume that k and k' are located 
near the K and K' points, respectively. When a uni- 
form magnetic field is applied, the path integration can 
be estimated, in terms of the magnetic length Is and the 
lattice constant a, as A{ri) ■ dri ^ o? For 



K ~K'\a, 



(A3) 



the intervalley mixing is negligible, since ($?|iJ|<I)^) ~ 

(t/N)Y,'e'^-^'^~'^'-^^ cx 6{k~k'). For the uniform field 
(perpendicular to the graphene sheet) of strength 10 T, 
one can find a/{l%\K-K'\) ~ 10"'' so that the condition 
(IA3p is achieved, which is why the intervalley scatter- 
ing can be ignored in current experimental studies. The 
intervalley mixing becomes important in a very strong 
magnetic field - lO'^ T, where a/{l%\K - K'\) - 0.1. 

In the same way as above, one can find the condition 
when the intervalley mixing is negligible in a nonuniform 
field with constant gradient A = VB(f). For the mixing 
to be ignored, the maximum value of the field must satisfy 
the condition (jASp . In addition, the gradient A should be 

much smaller than 10** T- A , since the path integration 
{e/h) Jf-, A ■ dri in Eq. (|A2[) becomes comparable to K ■ 

Ra-K'-Rb when a = Ao ~ 1.7x10'' T-A~\ In current 
experimental studies, the gradient is much smaller than 
Ao so that one can ignore the intervalley scattering. 



a! 



^ a - 



d 



Bo 



-Bo 



Bo 



^ c 



b? 



b° 



FIG. 6: The same as Fig. [3a) but with detailed labeling of 
edge channels (represented by solid arrows) at its four scat- 
tering points p £ {a,b,c,d} (dashed boxes). At each point p, 
there are two incoming channels with amplitude (pi,P2) a-nd 
two outgoing ones with (pfjPp). 



edges arid the magnetic edge states along the boundaries 
ab and cd. The scattering between the edge channels oc- 
curs at four scattering points p € {a,b, c, d}. Each point 
p has two incoming channels with amplitude {p\,P2) and 
two outgoing ones with {p'i,P2)', for example, the two 
incoming channels to the point a are one right-going and 
the other left-going channels along the upper ribbon edge, 
while the two outgoing channels from a are the two mag- 
netic edge channels along the line ab (see Fig. [6]). At each 
point p, we introduce a scattering matrix Sp which links 
the amplitudes of the incoming and outgoing states, 



P? 



M pI 



P\ 
P\ 



(Bl) 



The scattering matrix has a general form of 2 x 2 
unitary matrix, 



^11 ■^11 

'21 ''22 



1 - 



1 - t^e-^'^p 



(B2) 

On the other hand, while edge channels propagate from 
one scattering point p to its neighboring point p', they 
acquire phase accumulation ipppi. As a result, one has 



APPENDIX B: TRANSMISSION THROUGH A 
GRAPHENE RIBBON INTERFEROMETRY 

In Appendix [B] we derive the transmission probabil- 
ity through the inteferometry in Fig. [5Ka), based on the 
scattering matrix formalism. The resulting expression in 
Eq. (|B7p can describe the Aharonov-Bohm effect of the 
interferometry. One can easily obtain the transmission 
probability for other setups with different edge channel 
configurations by slightly modifying the derivation. 

The interferometry is in the integer quantum Hall 
regime so that its electron transport can be described 
by edge channels, such as the edge states along ribbon 



(4,4) 



e'^-(a?,aO), 



= (c?,c?). 



4 

02 



(B3) 
(B4) 
(B5) 
(B6) 



By combining the relations (jBip - (|B6|) and by putting 
o}^ — \ and (\ — 0, one can obtain the transmission 
probability T — \d^\^ of the edge state incoming from 
the source (the left of the inteferometry) to the drain 
(the right). 



T = 



(B7) 



8 



where = yl~t^, ip contains the Aharonov-Bohm 

phase as weh as the dynamical phase accumulated along 
one circulation of the closed loop abed, 

^ ^ ^ba+ Vcb + fdc + fad + 4''a- 4>b + 4>c - ^d 



+ 



p£ {a,b,c,d} 



(B8) 



a - 0a + 0& + 0'a - and /3 = <?!)c + (/^d + 0c - 4>'i- The 
transmission T can describe the Aharonov-Bohm oscilla- 
tion of the Dirac fermions in the setup of Fig.[Ha), as a 
function of B. Note that in general the scattering matrix 
Sp contains the information of the scattering between the 
K and K' valleys.-^ 

We close this appendix by analyzing Eq. jBTj for a 



simple case of ta = it = tc = td = t = 1/V2. In this 
case, the transmission probability can be simplified as 



T = 



cos^ # cos^ 4 



1 + 2 sin f sin f cos(v? - f - f ) + sin^ f sin^ f 

(B9) 

This result shows a usual form for the Aharonov-Bohm 
interference, except for the factors cosq;/2cos/3/2 and 
sin a/2 sin /3/2. From the factors, one can see that there 
appears no interference whenever a = tt or /3 = tt. It 
happens when destructive interference occurs during the 
propagation from one ribbon edge to the other through 
the two magnetic edge channels along ab or cd in such 
special cases. 
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